Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Gauge Theory for Quantum Spin Glasses 

Satoshi MORITA^, Yukiyasu Ozeki^ and Hidetoshi NiSHiMORi^ 

^Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro-ku, Tokyo 152-8551 
"^Department of Applied Physics and Chemistry, The University of Electro- Communications, 

Chofugaoka, Chofu-shi, Tokyo 182-8585 

The gauge theory for random spin systems is extended to quantum spin glasses to derive 
a number of exact and/or rigorous resuhs. The transverse Ising model and the quantum 
gauge glass are shown to be gauge invariant. For these models, an identity is proved that the 
expectation value of the gauge invariant operator in the ferromagnetic limit is equal to the 
one in the classical equilibrium state on the Nishimori line. As a result, a set of inequalities 
for the correlation function are proved, which restrict the location of the ordered phase. It 
is also proved that there is no long-range order in the two-dimensional quantum gauge glass 
in the ground state. The phase diagram for the quantum XY Mattis model is determined. 
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1. Introduction 

The problem of spin glasses has been attracting continued attention.^' ^ After the pioneer- 
ing work by Edwards and Anderson (EA),^ Sherrington and Kirkpatrick have exactly solved 
the mean-field model by assuming the replica symmetry.'^ Parisi has proposed the replica 
symmetry breaking solution and established the theoretical picture that the low-temperature 
phase is composed of infinitely many stable states with ultrametric structure.^ 

A topic of active investigations in recent years concerns the properties of short-range 
systems. Numerical studies have provided strong evidence for the existence of the spin glass 
(SG) transition for both the EA model^"^ and the XY gauge glass^'^"^^ in three dimensions 
but against it for the two-dimensional EA model. '^'^ For the two-dimensional gauge glass, 
although the long-range SG order has been denied rigorously, it is still possible that the 
system has a quasi long-range order in which the SG correlation decays in a power low. 
There remains the controversy about the existence of this order: some numerical studies have 
supported the absent of a finite-temperature transition^^' but some groups argue against 
such a conclusion. "'^'^'^^ 

Analytical calculations for spin glasses in finite dimensions are difficult because of ran- 
domness and frustration. However, a method using the gauge symmetry of the system is 
well-known as a powerful technique. ^^'■^'^ This method provides various rigorous results, for 
instance, the exact internal energy and an upper bound for the specific heat in the special 
region of the phase diagram. Another noteworthy result is a set of inequalities for the corre- 
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lation function, which restrict the topology of the phase diagram. In this relation, it has been 
suggested that the phase boundary between the ferromagnetic (FM) and SG phases is vertical 
by modifying the probability distribution.^^ These results are generalized to a wider class of 
systems including the usual Ising SG and the Zq and XY gauge glasses. 

Although the gauge theory provides us with surprising results, its targets have so far 
been limited to classical spin systems. In the present paper, we generalize this theory so that 
it applies to quantum spin systems. A difficulty of this generalization is the fact that we 
must define the gauge transformation for spins without violating the commutation rule. We 
circumvent this problem by using a rotational operator on the Hilbert space as the gauge 
transformation. 

This paper consists of six sections. In the next section, we formulate the gauge transfor- 
mation in two quantum spin glasses, the transverse Ising model and the quantum gauge glass 
(QGG), and show that these models have gauge symmetry. In §3, we prove an identity for a 
gauge invariant operator. This identity is valid even when the system parameters depend on 
time following the time-dependent Schrodinger equation. In §4, we derive a set of inequalities 
for correlation functions and order parameters. These results restrict the location of the FM 
phase or the Kosterlitz-Thouless (KT) phase in the phase diagram. In §5, we extend these 
inequalities to the ground state. The resulting inequalities for the order parameters show that 
the FM order does not exist at zero-temperature in the two-dimensional QGG. In §6, we 
consider the quantum XY Mattis model and determine its phase diagram. The last section is 
devoted to summary. 

2. Gauge Transformation for Random Quantum Spin Systems 

2.1 Transverse Ising model 

First, let us consider the random-bond Ising model in a transverse field. The Hamiltonian 
for this model is written as 



where af is the Pauli matrix at site i. Although we treat spin-1/2 systems in this paper, 
one can straightforwardly generalize all the results to spin-S" systems. There is no restriction 
in the spatial dimensionality or lattice structure. The exchange interaction Jij is a quenched 
random variable. One of the useful probability distributions for is the binary distribution 




(1) 



P{Jij) =p5{.hj - J) + (1 - p) 5{Jij + J). 



(2) 



It is convenient for later arguments to rewrite this distribution as 



(3) 
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where Tij = Jij/ J is the sign of the exchange interaction Jij. Another useful distribution is 
the Gaussian distribution 

nJ^.) = ^--p(-^-^^'^)^ (4) 



V2^ " V 2 J2 
where Jq and denote the average and variance, respectively. 

For quantum spin systems, the classical gauge transformation, which simultaneously 
changes the sign of all components, is not valid because the commutation rule [af, a^] = 2ia^ is 
changed to [o"f , a]'] = —2ia?. Thus we define a gauge transformation for spins using a unitary 
operator as 



U:af^ GafG-\ G = JJ ^i, Gi = < 



exp(-|:af) (^. = -1) 



(5) 



where is a classical gauge variable at site % and takes two values ±1. If = —1, a\'^ —af'^ 
and -^crf- Equivalently we can write 

U:iat,alat)^ia!,^iCTl^iat). (6) 

A difference of gauge transformations between quantum and classical systems is the transfor- 
mation rule of af. 

The gauge transformation for the bond variables {Jij} is the same as in classical systems, 
namely 

V : Jij ^ Jij^i^j. (7) 

The transverse-field term in eq. (1) does not change by the gauge transformation. 

The Hamiltonian (1) is clearly invariant under the successive operations of V and U: 
{UV)H = H. However, the distribution function of bond configuration is changed, for the ± J 
Ising model, as 

It is important for the following argument that this transformed distribution is proportional to 
the Boltzmann factor of a classical system. Similarly, the Gaussian distribution (4) is changed 
as 



Jfj + •^0 1 f Jo 



To simplify the arguments, we focus on the binary distribution (3) and (8), hereafter. It is 
straightforward to apply the same methods to the Gaussian distribution. 
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2.2 Quantum gauge glass 

Next, wc consider the quantum gauge glass (QGG). Similarly to the transverse Ising 
model, the gauge transformation is defined by the rotation operator. 

To properly define the quantum version of gauge glass, let us first consider the Hamiltonian 
of the classical gauge glass (CGG), 

Hc\ = -J^ cos(^j - - cjij). (10) 

This Hamiltonian can be rewritten using the spin vector composed of x and y components 
Si = (5) ^^'^ rotational matrix in the XY plane, R{e) = (^?^^ ~c^e^) as 

Ha = -jJ2'SiR{uij)Sj. (11) 

{ij) 

Thus this model can be quantized straightforwardly by replacing the elements of the above 
spin vectors by the Pauli matrices. The Hamiltonian of the QGG is therefore written explicitly 
as 

H = -jJ2{cosUi^ + - sina;,,- (af aj - afa]) } . (12) 

The phase factor ujij G [0, 2ir) is a quenched random variable whose probability distribution 
is of cosine type 

^Kp cos u>ij 

a periodic Gaussian (Villain) type 

^K) = Y^ 2^ exp( J (14) 

n=—oo ^ ' 

or a binary type 

P{LOij) =p5{L0ij) + {l-p)S{LOij -tt). (15) 

Equation (12) is a special case of the XY model with Dzyaloshinskii-Moriya interactions. 
This Hamiltonian is written as 

H = - j'^cTi- aj - '^ Jij{ai X crj)z, (16) 

where the second term is the random Dzyaloshinskii-Moriya interaction. If we set new param- 
eters, 

4 = y/j2 + j2, u;ij = -tan-^(^:^y (17) 
the above Hamiltonian is rewritten as 

H = -Y1 4- {cosu;., (af + afa]) - sinu^j [afa] - a\a^) } . (18) 
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This is equal to cq. (12) except that the interaction depends on bonds. 

In the CGG, the gauge transformation for spins is defined by the shift of spin variables as 
4^1 ^ 4'i~ i'i-, where V'j £ [0, 27r) is the gauge variable. Using the same notation as in eq. (11), 
this transformation is expressed as 

U:Si^ R{-^i)Si. (19) 

Thus we use this definition of gauge transformation for the QGG. Using a rotational operator 
on the Hilbert space, we define 

U:ai^GaiG-^ G = JJexp (-'-^aA. (20) 

The transformation rule for the transposed vector *(Ti is defined as 

U :'ai^'(TiR{i;i) = G'aiG-\ (21) 
The gauge transformation of random variables is the same as in the classical case, 

V : Uij Uij + il^j. (22) 
Under the gauge transformation UV, the Hamiltonian is invariant because 

{UV)H = ~jY,'^iR{^i)R{^i3 + ^3)Ri-^>3 = H, (23) 

(ij) 

where we used the property of rotation matrices, R{ip)R{(l)) = R((ip + (p). The probability 
distribution (13) is changed as 

This transformed distribution is proportional to the Boltzmann factor for the CGG. If we 
choose the Gaussian or binary type, the Boltzmann factor for the Villain or it J model appears, 
respectively. 

3. Identity for Gauge Invariant Operators 

The gauge symmetry of the Hamiltonian yields a useful identity for gauge invariant oper- 
ators. First, let us suppose that the system was initially in the perfect FM state \Fz) in the 
transverse Ising model. This state appears in the FM limit, T = 0, p = 1, h = 0. The gauge 
transformation operator G defined in eq. (5) changes this state as 

G|F,) = |0, I0 = l0i|02---I0iv. (25) 

If = +1, |0i denotes the state with up spin in the z direction, and if = —1, the spin at 
site i is down. 

Using the property of \Fz) in (25), we prove the following identity for a gauge- invariant 
operator Q which satisfies Q = {UV)Q or equivalently VQ = G~^QG, 

[iQ)F.] = [(Q)pel(i^.)] ' (26) 
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where (•)^ i{Kp) ^^e expectation value for the classical equilibrium state on the Nishimori 
line (NL), that is, 

(Q)p.(K,) = IVpe.(i^,)Q, Pc.(K,) = ^£Z!^. ■ (27) 

To prove the identity (26), we apply the gauge transformation for the configuration of 
randomness of cq. (7) appearing on left-hand side of cq. (26). This operation does not change 
its value because the transformation V of cq. (7) only changes the order of the summation 
over Tij. Thus the left-hand side of eq. (26) is rewritten as 

where we used the assumption that the operator Q is gauge invariant, VQ = G~^QG. Since 
the expectation value on the left-hand side does not depend on ^, the summation over ^ and 
division by 2^ yield 

The last part of the right-hand side is rewritten in terms of Pci(-Kp) as 

J2 e^- S {(\Q\0 = j J] e^- S I (TVpei(i^p)Q) • (30) 



Therefore, we obtain 



cosh hp J 



Since Trpci(-fi'p)Q is invariant under the transformation V, this is identical to the right-hand 
side of eq. (26) 

The above result can be generalized to the case that the transverse field h{t) depends on 
time, following the classical example.^^ We consider the zero-temperature time evolution fol- 
lowing the Schrodinger equation. Using the time ordered product, the time evolution operator 
is written as 



(^-i J\{t')dt'^ . (32) 



?7t = rexp (^-i 

Since the time dependence of the Hamiltonian does not invalidate the gauge symmetry, this 
operator is also gauge invariant 

{UV)Ut = G{VUt)G-'' = Ut. (33) 

Examples of the gauge invariant operator include the transverse magnetization af{t) = 
u}afUt, the autocorrelation function a^{0)a^{t) and the interaction term Ho{t) of the Hamil- 
tonian (the first term on the right-hand side of eq. (1)). 
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For the QGG, we can prove a similar identity 

[iQ)Fj = [{Q)pAk,)] ■ (34) 

Here one should note that Pc\{Kp) is different from the normal density operator. If we choose 
the cosine-type distribution (13), Pc\{Kp) is defined by the Boltzmann factor for the CGG as 

Tl^e^^ cos(a;i3-V'i+V'j) I 



PciiKp) 



(35) 



where = G\Fx) and Tr^ stands for integration over ijji from to 27r. Since the state vector 
l-^) does not diagonalize the Hamiltonian for the quantum gauge glass, Pci{Kp) is not equal 
to the density operator 



P{Kp) = 



(36) 



The identities (26) and (34) show that the expectation value of gauge invariant operator in 
the FM limit is equal to the one in the classical equilibrium state on the NL. The equivalence of 
the two states has already been proved in the dynamical gauge theory for classical systems. ^^'^^ 
The present results are generalization of these dynamical cases to quantum systems. Note that 
the zero-temperature time evolution for quantum systems is deterministic in contrast to the 
stochastic dynamics for the classical SG. 

4. Correlation Function and Order Parameter 

Using the identities proved in the previous section, we can derive a class of inequalities 
for the correlation function. First, we treat the transverse Ising model. Since the correlation 
function is not invariant under the gauge transformation, let us consider the following gauge- 
invariant quantity 

Q = ala]{a!a^^)K,h, (37) 

where {■)K,h denotes thermal average with temperature = J/K under a transverse field 
h. Substitution of the above quantity into eq. (26) yields 

Here (CiCj)iCj, is the correlation function for the classical Ising system with the same configu- 
ration {Tij} and no external field. By taking the absolute value of both sides of this equation, 
we find 



< 



(39) 



where we used the fact that correlation function ( crfaf ) does not exceed unity. Similarly. 



K,h 



we can prove 



< 



(40) 



7/15 




Fig. 1. The phase diagram of the transverse Ising model. The paramagnetic (PM), the ferromagnetic 
(FM) and spin glass (SG) phases meet at the multicritical point (MCP) in the plane h = 0. The 
dashed line at p = pc sets a bound for the existence of the FM phase also for h ^0. 
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If the probability p of the FM interaction is less than the critical probability pc at the 
multicritical point for the classical Ising system, the right-hand side of the inequality (39) 
vanishes in the limit |i — j| — > oo. Thus the correlation function for the transverse Ising model 
on the left-hand side is also equal to zero. Therefore the region of the FM phase is restricted 
to the range p > Pc (Fig. 1). Since the transverse field represents quantum fluctuations, the 
correlation function for the transverse Ising model should be reduced from classical system 
with h = 0, which is the physical origin of the above-mentioned restriction. 

Next, let us consider how to define a correlation function of the QGG which has convenient 
properties for the gauge theory. In the CGG, a useful correlation function is defined in terms 
of an exponential function as e'*^'^*"''^^^. This is rewritten using the notation of eq. (12) as 

^KcPi-cPj) = Si-Sj-i {Si X Sj), = ^SiSj - i*SiR{-iT/2)Sj. (42) 
This motivates us to define a correlation operator for the quantum gauge glass as, 

li3 = + ^t<^]) - i(^f - ^^^J)- (43) 
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Ujij = GjijG ^ = e" 



(44) 



llj — <J /jj»_r — c - 7ij. 

The first factor is the same as the correlation function appearing in the classical gauge glass 
except for a minus sign. 

Usually, the correlation function for XF-like systems is defined as the expectation value of 
afaj+afa^, which corresponds to cos((?!)i — 0^) for classical systems. However, this expression 
is not useful for the gauge theory because it does not separate into gauge variables and spin 
operators after transformation by U as 



In addition, one can easily prove that 

since the QGG is invariant under the following transformation: 

Thus we choose to discuss the correlation operator (43). 

To derive an inequality similar to eq. (39), it is useful to consider 

Q = 7ij{lij)K- 



(45) 

(46) 
(47) 



(48) 



It is easy to prove that this quantity is gauge invariant because of the property of given 
in eq. (44) . Substituting this quantity into eq. (34) , we obtain 

cl 



.-i(V'i-V'j) 



(49) 



where (•)^^ stands for the expectation value for the CGG on the NL. By taking the absolute 
value of both sides of this equation, we find 

\[{li3)K\\< 



i{ipi-'4>j)Y 



I 



Similarly, we obtain 



[sgn((7ij)if )] I < 



cl 



K 



K 



> 1. 



(50) 



(51) 



(52) 



Since the lower critical dimension di is two for continuous spin systems, for d > 2, we 
expect a FM phase to exist at low temperature under small randomness. In this case, a 
two-spin correlation function tends to the square of magnetization when the two spins are 
sufficiently separated, 

[(7y ] rn{K, Kpf, \i - j\ 00. (53) 
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The right-hand side of inequahty (50) is estimated as follows: 



cl 



< 



-iiipi-itj) 



Kr. 



(K-Jl 



OO 



(54) 



= q''\Kj„Kpf=rrf\Kp,Kpf, 

where we used the identity rrf^ = q'^^ resulting from the gauge theory on the NL for the CGG. 
Therefore we obtain 

m{K, Kpf < rrf\Kp, Kp). (55) 

If the parameter Kp is smaller than the critical point Kp^ for the classical system, the right- 
hand side vanishes. Consequently, the FM phase for the quantum gauge glass lies in the region 
satisfying Kp > K^K This result is consistent with the intuitive picture that quantum effects 
reduce long-range order. 

If the spatial dimensionality of the system is equal to the lower critical dimension, d = 2, 
there is no long-range order but quasi long-range order. The Kosterlitz-Thouless (KT) phase^^ 
exists when both K and Kp are sufficiently large. The ordering tendency of the KT phase is 
observed by the correlation length ^ in the paramagnetic (PM) phase as 



K. 



~ e 



-\i-j\/U{K,Kj,) 



Kr. 



cl 



\i-3\/£.'rL[K^,Kj,) 



(56) 

(57) 
(58) 



From the square of the inequality (50), the limit |i — j| — > oo yields 

U{K,Kp)<2it{Kp,Kp). 
Thus, the KT phase for the QGG is also restricted to the region Kp > KpK 

5. Ground State Property of Quantum Gauge Glass 

Next we consider the ground state property of the QGG. It is necessary to consider the 
transformation rule for eigenstates of the Hamiltonian. Let us denote an eigenstate of the 
Hamiltonian with the eigenvalue x by |x;u;), 



H\x;u}) = X \x;u)). 
Since the invariance of the Hamiltonian (23) can be rewritten as 

G{VH)G-^ = H, 

we have an eigen-value equation 

(VH) G-^\x;u;) =x G-^\x;u;). 



(59) 



(60) 



(61) 
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Note that the effect of the operator V is restricted to the inside of the brackets (•••). The 
state G~^\x;u) is an eigenstate of the gauge-transformed Hamiltonian (VH). Therefore, one 
can derive the transformation rule of the eigenstate 

V\x;u}) = G-^\x;uj). (62) 

This rule is derived straightforwardly, if the state is not degenerate. When the state is degen- 
erate, the rule is not unique. However, it is reasonable and has no problem, if we use this rule 
as the transformation rule. 

Now, we consider the ground sate. Let us denote the ground state of the Hamiltonian by 
|g;a;). The above transformation rule leads to the gauge transformation of the average of any 
operator Q in the ground state; 

V{g;u;\Q\g;u) = {g;u;\GiVQ)G-'\g;u:) . (63) 

The thermal average {■ ■ ■)k in equations derived in the previous section can be replaced by 
the ground state expectation value. For example, 

cl 



[(g;a;|7ij|g;a;)] 



(64) 



is derived instead of eq. (49), which provides the inequality for the order parameters in the 
ground state, 

m(oo, Kpf < rrf^{Kp, Kp) (65) 

instead of eq. (55), and 

U{^,Kp)<2C{Kp,Kp) (66) 

instead of eq. (58). 

In two dimensions, it has been shown that the FM long range order exists in the ground 
state of the pure quantum XY model {Kp = +00).^^ However, in the disordered regime 
{Kp < +00), the FM order must disappear since the FM order in the CGG model, the right- 
hand side of the inequality (65), does not exist. The only possibility in this regime is the 
existence of the KT phase, which is consistent with the inequality (66). 

6. Phase Diagram of Quantum Mattis Model 

In this section, wc introduce and discTiss the properties of the quantum XY Mattis model 
which has no frustration. Using the gauge transformation, one can obtain the phase diagram 
for this model explicitly. One of the phase boundaries is determined by the critical point of 
the pure quantum system and the other by that of the classical one. This is an important 
difference from the classical non-frustrated systems. 

Let us locate a quenched random variable at each site and define the phase factor ujij 
as LOij = Uj — (jJi- The Hamiltonian for the quantum Mattis model is defined in terms of the 
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pure quantum XY model, Hq, that is, 

H = G^HoGz' , G^ = ll exp (- ^'^f ) • (67) 

Thus, the ground-state energy is always the same as that of and the ground state is obtained 
by operating G^^ on the ground state of the pure system. This system has no frustration in 
this sense. It is easy to show that the Hamiltonian for this model is invariant under the gauge 
transformation, (UV)H = H, where the gauge transformation for the configuration is defined 
as tVi ^ uji — ipi. 

Using eqs. (44) and (67), we immediately obtain 

where the angular brackets on the right-hand side denote the thermal average with respect 
to Hq. Note that the correlation function for the pure system {jij)^ ^ does not depend the 
quenched variable coi. Here, we assume that the distribution function for the quenched random 
variable is proportional to the Boltzmann factor of the pure classical XY model. 

Hereafter, the subscript and the superscript cl stand for pure and classical systems, respec- 
tively. Then the configuration average is equal to the thermal average for the pure classical 
XY model with coupling Kp. Consequently, we find 

[(7.,);,]=(e^("^-"^-))^'^^(7..)o,^. (70) 
Similarly, from eq. (68), the spin-glass correlation function satisfies 

21 2 
K7u)xl J = {li3)Q,K ■ (71) 

Taking the limit \i — j| ^ oo of eqs. (70) and (71), we obtain 

m{K, Kp) = ml\Kp)moiK), (72) 

q{K,Kp) = mo{Kf. (73) 

For d > 2, there are three phases, PM, FM and Mattis spin-glass (MSG) phases. Figure 2(a) 
shows the phase diagram. The phase boundary between PM and other phases is at K = Kqc 
and the one between FM and MSG is at K = K^l- 

li d = 2, the correlation length determines the phase structure. Prom eqs. (70) and (71), 
we find 

' +^ (74) 



U{K,Kp) ^^\Kp) UK) 

UK,Kp) = ^-^, (75) 
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Fig. 2. (a) The phase diagram for the quantum XY Mattis model in d > 2. There are three kinds of 
phases, the PM, the FM and the SG. (b) The phase diagram in d = 2. There exist three kinds of 
phases, the PM, the uniform KT (UKT) and the random KT (RKT). 



where ^q{K,Kp) denotes the spin-glass correlation length. Thus, similarly to the d > 2 case, 
three phases exist: (i) < oo and < oo: paramagnetic phase (PM), (ii) = and 
= oo: uniform KT phase (UKT) and (iii) < oo and = oo: random KT phase (RKT). 
The phase diagram is shown in Fig. 2(b). 

We note again that the location of the horizontal phase boundary is determined by the 
critical point of the quantum pure system, -fCoo and the vertical one comes from that of the 
classical pure system, K^^. 

7. Conclusions 

In this paper, we have investigated quantum spin glasses using the gauge theory. First, 
we considered the transverse Ising model and the QGG. To construct the gauge theory, we 
defined the gauge transformation by rotational operator on the Hilbert space. It is essential 
that interaction of these models is written in term of one or two components of spin operators. 
If a system has Heisenberg-type interactions, we can not define a gauge transformation which 
satisfies the commutation rule. 

Using the gauge theory, we obtained mainly two results. One is the identity for gauge 
invariant operators. The FM limit state and the classical equilibrium state on the NL provide 
the same expectation value for gauge invariant operators. We note that this result remains 
valid when we introduce time evolution following the Schrodinger equation. This result has 
already been pointed out for classical spin glasses with stochastic dynamics.^^'^^ We have 
shown that the same applies to quantum spin systems. 
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The other result is a set of inequahties for the correlation function. These inequalities show 
that the correlation function for the quantum model never exceeds the classical counterpart 
on the NL. The corresponding classical system is determined by a transformation rule for 
probability distribution. As a result, the order parameter is smaller than the square of the 
classical one. Therefore the FM phase (or the KT phase) should lie within the corresponding 
classical one. This is natural intuitively since quantum effects reduce ordering tendency, but 
to prove it rigorously is a different and quite a non-trivial problem. Moreover, these results 
are valid even if the system is in the ground state. Thus FM long range order vanishes in the 
two-dimensional QGG although the ground state of the pure quantum XY model has FM 
order. 

Next, we determined the phase diagram for the quantum XY Mattis model. This model is 
not a real SG because of lack of frustration. It is interesting, nevertheless, that both quantum 
and classical phase transitions occur in a single system, which may serve as a starting point 
for investigations of more realistic quantum spin glasses with frustration. 
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